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Abstract 

Confinement is explained via field correlators, and the latter are 
calculated via gluelumps. Behavior of gluelump Green's function at 
small and large distances yields gluonic condensate and vacuum cor- 
relation length respectively and allows to check the consistency of the 
whole picture. 

1 Introduction 

Field theory based on Abelian and nonabelian gauge fields is now the foun- 
dation of physics of strong and electroweak interaction. Here the gauge co- 
variance in general and gauge invariance of physical amplitudes is the corner 
stone and the most general principle, which allows to fix the form of equa- 
tions and resulting amplitudes even before any dynamical input. In the case 
of QCD this fact is even more important, as was stressed in [H [2]. Below we 
show, that confinement makes gauge invariance and all relations, connected 
to it, the most important property sine qua non, and in turn, gauge invari- 
ance principle allows to build theory of confinement in good agreement with 
all existing data. 

It will be shown below, that gauge invariance requires each particle trav- 
elling from point x to point y to develop its own gauge factor $c(a;, y) (the 
Schwinger line, or parallel transporter) pertaining to the path C{x,y). For 
neutral (color-singlet) systems the set of paths {Ci{x, y)} forms a closed loop. 
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a simple Wilson loop W{C) for two opposite charges, or more complicated 
closed loop with A^^ branches for the Nc charges in the group SU{Nc). 

As will be shown (see also [3] and |1] for the review) all dynamics in the 
neutral system is given by the vacuum-averaged Wilson loop {W{C)), and 
since all physical amplitudes correspond to the gauge invariant and hence 
neutral systems, therefore {W{C)) contains all the clues to confinement and 
in general to the dynamics of strong interaction. 

Here one can distinguish confining theories, like QCD, and nonconfining 
ones, like electroweak (EW) theory. In the first case gauge invariance and 
confinement make it impossible to consider only one part of color neutral 
system, e.g. propagator of quark or gluon. Formally, gauge invariant vac- 
uum average of this colored object is zero, physically, the confining string, 
which appears in {W{C)) and connects this object to its neutralizing partner, 
strongly governs its motion and cannot be disregarded. Therefore all models, 
based on vacuum averages of gluon or quark propagators in any gauge, are 
irrelevant to the nature. 

In nonconfining theories the situation is different. One can compose a 
system of an electron here and proton on the moon and consider the motion of 
electron in the gauge invariant language of {W{C)), but the latter factorizes 
into a product of propagators, however some initial and final conditions of 
electron should be fixed to avoid divergencies. 

A similar situation occurs in the deconfined phase of QCD [S] (see [B] for 
a review), where (apart from magnetic confinement) {W{C)) also factorizes 
and dynamics is given by Wilson lines rather than Wilson loops. 

Thus the basic problem is first the analysis of {W{C)) and, second, self- 
consistent calculation of {W{C)) or its ingredients. The first part is done 
using the gauge invariant cumulant expansion of {W{C)) in terms of field 
correlators. As was shown in |S1 IH E], for QCD the lowest term (quadratic) 
{F{x)F{y)) saturates {W{C)) with few percent accuracy, which is supported 
by the Casimir scaling measurements on the lattice both in the confined [8] 
and deconfined phase [9]. The same Casimir scaling helps to put a strong 
upper limit on the presence of topological charges or adjoint fluxes [7j as a 
possible source of confinement0. At the same time, Casimir scaling is well 
supported by the Gaussian quadratic term {F{x)F{y)) and the latter yields 

^ Note, that magnetic monopoles in QCD contribute largely to higher cumulants and 
hence violate Casimir scaling (see [TOt) while in the weakly-coupled dilute 3d SU{N) 
Georgi-Glashow model the Casimir scaling holds to a good accuracy [11], implying that 
higher cumulants are suppressed at low density. 
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a beautiful picture of linear confinement in good agreement with all lattice 
measurements and hadron physics. 

Here comes the second part of problem: which field configurations are 
behind the Gaussian term, and can one calculate these correlators selfcon- 
sistently, so that all QCD can be derived from the one constant, fixing the 
scale of our world. We show below, following [12], that this is in principle 
possible, and one can reexpress some physical observables, like string ten- 
sion cr, through others like Aqcd (however still with low accuracy). More 
importantly, one can establish the objects, which are behind confinement 
and characterize full vacuum structure: those are known as gluelump^, and 
consist of one or two gluons propagating in the field of adjoint static source. 
Mass of these objects was calculated analytically before [U] and its inverse 
gives the correlation length of the vacuum, while the short distance behav- 
ior yields gluonic condensate, and finally, the integral of gluelump Green's 
function yields string tension. Thus confinement is connected to gluelumps, 
and all other vacuum properties can be calculated via gluelumps. Moreover, 
the small correlation length A (due to large gluelump mass M, A ~ ^/M) 
explains, why the Casimir scaling holds in QCD - the cumulant expansion is 
actually a series in powers of the dimensionless parameter cxA^ = 0.05 , and 
the first (quadratic) term is dominant. Below we give a short exposition of 
Wilson loop and cumulant expansion in section 2, gluelumps are introduced 
in section 3, and exploited in section 4 to display consistency of the whole 
picture. Summary and conclusions are given in section 5. 



2 Wilson loop and field Correlators 

The Green's function (propagator ) of quark can be written in the Fock- 
Feynman-Schwinger Representation (FFSR) (see [16] for a review) as the 
quantum mechanical path integral (we consider Euclidean space-time through- 
out the paper) 

S{x,y) = {m+Oy^ = (m-D) / ds{Dz)^ye~^ PAewiw / At^dz^,)p^{x,y, s), 

Jo Jy 

(1) 

^The name and the first calculations in SU{2) lattice theory are due to I.H.Jorysz and 
C.Michael [13]. The detailed SU{3) calculations have been done in [M] 
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where K = ms + \ dr and = PA^y^piig !y ^iJidz^j)p(j is the 

Schwinger hne (parallel transporter) with spin insertions, and 



p„{x,y,s) = Ppexp 



9 [ a^yF^^{z{T))dT , (2) 



The last factor contains spin operators 

^M-^M- = ( cr . E cr ■ B ) ' 

and will be important for calculation of gluonic condensate in section 4; it 
also yields all spin-dependent interactions in hadrons, see [17] for a review 
and references. 

In an analogous way one can construct FFSR for a gluon in the external 
and vacuum fields 

G{x,y)= r ds{Dz),ye-^PAe'^y^''''p^ (4) 
Jo 

where implies in the adjoint representation and the spin factor can 
be written either in the form ([3]) with the gluon spin operator S instead of 
cr, or simpler, as 

PS = [PFexp(2z(7 r Fx^{z{T))dT)]^,. (5) 
Jo 

We omit in this section the spin parts of quark and gluon Green's function, 
since they do not define the confinement picture at this stage. 

For a white hadron, consisting of qq or two gluons one can write the total 
Green's function as 

iC)W{C,y) (6) 

and the Wilson loop operator for the closed loop C^y is 

W{C) = {PAexpitg f_ A^dz^)), (7) 

where the brackets imply vacuum averaging. The use of the nonabelian 
Stokes theorem and cluster (cumulant) expansion fT8] yields 



=tr,exp|f;M! J dsil) J rf.(n)((F(l)...FH))|, (8) 



4 



where ds{k) = ds^^^y^ is the surface element, and 

F{k) = F^^y^{uk,Xo) = ^{xo,Uk)F^^y^{uk)^{uk,Xo) (9) 

while y) is the parallel transporter. Field correlators {FC){{F {!)... F{n))) 
depend in general on the chosen point xq, which is convenient to choose 
somewhere on the minimal surface bounded by C. The great simplification 
for QCD is that FC decrease so fast with distance, as we shall see later, 
that xo)$(xo, Mj) in f l8|9|) can be replaced by straight line transporter 

Uj), and moreover, the whole series is fast converging as (crA^)", where 
A is the correlation length, {F{x)^F{y)^) ~ exp (~ ^^a^^ ) • Therefore one can 
keep only the Gaussian correlator and moreover, since A is much 

smaller than all other characteristic quantities, e.g. hadron sizes, inverse 
masses of lowest hadrons etc., one can say that the QCD vacuum configura- 
tions correspond to the Gaussian white noise. 

As a result one can use only the lowest correlator 

D,,,M = {trj^F^,{x)^{x,y)Fx.{ymy,x)), (10) 
which can be expressed in terms of two scalar functions D{z) and Di{z) 



OZf, 



Di{z). (11) 



Note, that (fTOj) . (ITT]) is written for fundamental charges, for charge of 

c{j) 

C(fund) ' 



irrepr. j one should multiply by the ratio twP^-, where C(j) is the quadratic 



Casimir operator. 

Now all the spin-independent dynamical properties can be calculated from 
static potentials, the confining Lorentz scalar V{r), 

= ^^I"- - A)<iAf d.i)(VA^) (12) 

and the Lorentz vector Vi(r) 

Viir) = \d\ rfz/Z}i(VA2 + z/2) (13) 
JO Jo 
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and the string tension (and confinement) is due to D{z): 



Spin-dependent potentials can be found as well for quarks of any mass or 
gluons without l/Mg expansion, keeping in ^ also spin terms, see [I7] 
for a review. The interaction (fT2i) . ( |T3|) was used in the framework of the 
relativistic string Hamiltonian [19] and spectra of mesons [20], glueballs, [2T] . 
hybrids [22] have been computed in excellent agreement with lattice and 
experiment with minimal, input (current quark masses, a and as without 
any additional parameters) see [23] for reviews. 

At this point one could wonder how to check the notion of Gaussian 
vacuum? The answer is very simple: Gaussian correlators D and Di are 
proportional to [gFy, i.e. to the charge squared, and hence for higher color 
representations j they are proportional to the Casimir operator C2{j), known 
for the group SU{N,), C2(fund) = C2{adj) = N, etc.. 

This constitutes the Casimir scaling, and the numerical check was done 
on the lattice [H |9] and compared with FC in [^. results of the lattice 
measurements in the confined and deconfined phases of QCD are shown in 
Fig. 1 (taken from [8]) and Fig. 2 respectively. One can see an exellent 
agreement of Ygq^j) = {—j:\nWj{C)) computed on the lattice for different 
j = 3,6,8,... with the Casimir scaled potential ^^^Vqq{3). (Note in Fig.l 
V and r are units of r^^ and tq respectively, vq ^ 0.5 fm.) 

In the deconfined phase only Vi{r,T) is nonvanishing (apart from col- 
ormagnetic correlators, see p] ) and it enters the Polyakov lines, Lj{T) = 
exp (~§|3y ^^^27''^'' ) • Oils can see in Fig.2 the lattice measured values 
-^fund(^), -^^adj(7') as fuuctious of temperature (points) vs our analytic calcu- 
lations with Casimir factors. 

Agreement is impressive. (Note, that at T < in confinement phase 
Lfund = 0, formally since one should replace Vi — )■ V^i(oo,T) -|- V{oo,T) and 
the latter is infinite). 



3 Gluelumps and Field Correlators 

To express FC through gluon propagators, one should consider gluons in 
Di_iu,XcT, i.e. gluons emitted in F^^{x) and absorbed in Fx^iy) and moving in 
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Figure 1: (from [8]). The potentials Vd in units of I/tq for all measured 
representations D = 3,8, 6, obtained at (3 = 6.2(ro = 0.5 fm, a = ro/6.1). 
The line for a representation D is the fundamental potential, multiplied by 
the ratio 



C2(D) 



C2(fund) • 



the background of vacuum fields (consisting of gluon fields again and with 
the source in $(x,?/)). This is the basic idea of mean field approach. It is 
important, that we can express the interaction of propagating gluons with 
background field again in terms of FC, and thus one obtains a closed scheme, 
consistency of it will be discusses in the next section. Expanding F^i, into 
Abelian and nonabelian parts, F^i, = {d^Ai, — d^A^) — ig[Afj_,A^] one can 
write 

D,.,xAx, y) = D^'l. + d^^Im + ^Sa. (15) 
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Figure 2: Shown on the figure are curves of Ladj (lower curve for all T/Tc) 
and Lfund (upper curve) compared to the lattice measured from [9j. In the 
T < Tc region the M{ds = 0.195) = 0.982 GeV gluelump mass was used. In 
the deconfinement region the fit was used with Tc = 270 MeV for Lfund and 
the Casimir scaled value for Ladj- 

where the superscript 0,1,2 refers to the power of g, coming from the term 
ig[Af,,A^]. 

Here D^^^^^ is connected to the one-gluon gluelump Greens' function 

DZ<^i^,y) = ^ |^|-G(^/)(x,2;) +perm| + a(°),,, (16) 

where the Ig gluelump Green's function is 

G%^\x,y) = {TTaA,{xMx,y)A,{y)). (17) 

and Tr^ implies summation over adjoint indices. 

^^i^uXa is of basic importance, since it ensures confinement via D{z) and 
is expressed via two-gluon gluelump Green's function G^'^^\z) 
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D^'i^,^{x,y) = {6,,6.,-6,J.x)D{x-y); D{z) = ^'^^^ ^^ G^'^\z), (18) 

where both one- and two-gluon gluelump Green's functions can be written in 
terms of path integrals [16] and finally expressed via eigenvalues and eigen- 
functions of relativistic string Hamiltonian \W\ . For gluelumps this Hamilto- 
nian was studied in [12] and results for lowest gluelump masses in comparison 
with lattice data [H] are (expressed via fundamental string tension cr/ = 0.18 
GeV2). 

jVf(i5) = 1.49 GeV, M(^5)(iat) < 1.7 GeV, 

M(29) = 2.61 GeV, M(29)(lat) = 2.7 GeV, (19) 

where M*^*^^ is the mass of gluelump with z = 1, 2 gluons. The same Hamil- 
tonian technic allows to find the asymptotics of Green's functions [12] 



D,iz) = ^WW|tf!l^^-Ma)|.| 1^1^(1,) ^ ^ (20) 



D{z) = ^'^^' 1) 0.1^2 -M(^^)kl M(25)|;^|»l (21) 
2 ■' 

Here is the value of the strong coupling constant at large distance. 
Correspondingly, one obtains the correlation length of D{z) and Di{z) 

^ = = 15 fa- (22) 

This should be compared with A, Ai calculated on the lattice previously 
[241, A ~ Ai ~ 0.2 fm, while a recent measurement [25], analyzed in [17] 
yields A ~ 0.1 fm, and the numerical HP^ formalism [26] yields A = 0.13 fm, 
Ai ^0.145 fm. 

Finally, a remarkable agreement with the values fl22p has been found in 
Ref. [27], where the full calculation of the path integrals representing Gj^^f-* 
and G^^^^ has been performed. 

Thus all data support the conclusion, that the correlation length is very 
small A, Ai ~ 0.1 fm and the QCD vacuum is therefore highly stochastic. 
Moreover, we have calculated FC in terms of nonperturbative parameters, 
like a, which can be calculated, as in f lT^ via FC themselves. This calls for 
the consistency check, to be done in the next section. 
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4 Consistency of the gluelump mechanism 



We start with small distances and with Di{z). 

The Ig gluelump Green's function can be written in FFSR [T6 



G';,^\x,y) = Tia / ds{DzUexp{-K){W;^,{C.y)), (23) 

<^ 

where the spin-dependent Wilson-loop with insertion of operators F^j, is 

KuiC.y) = PPf J A^dzx) exp f| (24) 

where exp F is as in 

With zero background {A^ — )■ 0) this is easily calculated and yields the 
lowest order perturbative result, Di{z) = then from (fT3l) Vi{r) = — 

The first nonperturbative (np) contribution comes, when one accounts 
for the small Wilson loop behaviour (second paper in [3]) 

W^;. = exp(-^G252), G2 = ^((F;,(0))2), 5 area. (25) 

As was shown in [12], insertion of (125 p in (125]) yields finally 



A(.) = l^^i±^MM + ^G,. (26) 

Note two important facts; 1) perturbative and np contributions in (126|) enter 
additively in this lowest order, this supports the original idea in |28j; 2) the 
sign of np part is positive and this contributes to the np shift down of the 
vacuum energy density. 

We turn now to the function D{z), connected via (fTSl) to the 2g gluelump 
Green's function. Here the small distance region is more complicated. First 
of all, the same small-area term fl25|) yields for D{z) a negative result! 

Diz)=A^Diz)+A2Diz), 

AM^) ^ (27) 

However the same spin-dependent Wilson loop (now with two gluons and one 
parallel transporter (Schwinger line)) in flMj) contains paramagnetic contri- 
bution {W^{1)W^{2)) ~ exp(-c(FF)), which has positive sign (see Fig. 
3). 
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Figure 3: Interaction of two gluon trajectories, depicted by broken lines, via 
correlator of operators F^^. 



A.Diz) = 2Nyhiz), h{z) ^ ^ In^ j (28) 

Here Aq is the point, which is above the asymptotic free region of as{X{z)), 
Ao ^ 0.24 fm (see Appendix of the last paper in [I2]). Note an interesting 
feature of fl28l) : for vanishing z, g^^z) ~ i^ij^ ^^"^ ^2D{z) tends to a positive 
constant, so that the total contribution, 

D{z -> 0) = A,D{z) + A,D{z) DiO) = ^DiXo) (^] . (29) 

2vr^ \Po J 

For Nc = 3 this yields -D(O) ~ 0.15D{Xq), and implies the humpbacked 
form of D{z), see Fig. 4, since it decreases exponentially at large z, see 

Note, that this exponential behavior D{z) = Da-{0)e~^^^^'^, if continued 
to zero, yields string tension from Eq. (fT^ a = 7rA^Do-(0), and -Do-(O) = 0.35 
GeVl Instead D{0) = fgGa and Ga from [28] (and this agrees with transition 
temperature Tc, as shown in [29]), G2 ~ 0.01 GeV^, which yields -D(O) ~ 0.01 
GeV^. Thus we understand qualitatively, that this discrepancy between -D(O) 
and -Do.(0) is due to the AF behavior of (125]) . illustrated by Fig.4. 

Now let us discuss the selfconsistency at large distances. We start with 
the asymptotic at \z\ > A, D{z) =t D^(0) exp(-|0|/A), A = 1/M^^^ with 
Drj{0) given in fl2T|) . and assume, that this form can be extended to z = 0, 
at least in the integral cr = ^ J (PzD{z), where small z are less important. 
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Figure 4: Qualitative behaviour of the correlator D{x) versus x, as predicted 
by (|20|) and fl27j) . fl28|) at large and small distances x respectively. 



As a result of this insertion of D{z) into ( fT4l) . one obtains an interesting 
relation 

0.l-87T'aliN'^-l)aj = ^ (30) 

where in = as(/i), and the scale corresponds to the average momentum 
(^ inverse radius) of the two-gluon gluelump. One can take /^o ~ 1 GeV, 
then fromf lSU]) one get^ 

An 

a,(^^l GeV) >0.4; «,(^) = (31) 

/3oln ^'^ + 



1? 

^^QCD 

where we have used the IR saturated form of a^, suggested earlier in [30], and 
derived for QCD with np background in [31] and tested in hadronic physics 
in [32] with the so-called background mass parameter Mb = 1 GeV. Eq. 
f l3T]) is actually the "dimensional transmutation relation", since it connects 
cr = 0.18 GeV^, as a basic scale parameter of QCD, yielding all masses and 
A, to as{fi), i.e. Aqcd- One gets from f l3T|) for three flavors n/ 

A(% = 3) > 0.29 GeV. (32) 



■^Note the sign > in (1311) due to the fact, that the extension of exponential form of D{z) 
to smaU z, overestimates the integral, cf. Eq (|29l) 
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One should take into account, however, that we have used Eq. fl31l) in the 
x-space, where A should be approximately 1.3 times larger, than Aqcd('^/ = 
3) ^ 0.27 GeV. 

Thus this large distance consistency check is reasonably satisfied. 
We end this section with two remarks: 

1) Till now only np contribution to D{z) have been discussed. Naively 

4 

the 2g gluelump Green's function yields perturbative result for D{z) ~ ^4, 
however, as shown in |33], all perturbative contributions to D{z) are exactly 
cancelled by those from higher correlators. 

2) D{z) and Di{z) have been obtained here in the leading approximation, 
when gluelumps of minimal number of gluons contribute: 2 for D{z) and 1 
for Di{z). In the higher orders of 0{as) one has an expansion of the type 

D{z) = D^^''\z) + cialD^^''\z) + C2alD^^''\z) + ... 

D^{z) = D^^''\z) + d^alD^^''\z) + ... (33) 

Hence the asymptotic behavior for D{z) will contain exponent of Mo^^|2;| too, 
but with a small preexponent coefficient. 

5 Summary and discussion 

We have shown, that every physical amplitude can be expressed through the 
corresponding Wilson loop, integrated as a path integral with the known 
weight. Therefore all perturbative and nonperturbative expressions for any 
amplitude can be derived from the corresponding expansions of the Wilson 
loop, e.g. the double-logarithm asymptotics of Sudakov type also follows from 
this representation [34j . On the other hand, Wilson loop can be exactly given 
by the cumulant series, containing all FC, and this is exact representation 
both for perturbative and np contributions. 

We have also demonstrated, that the np part of the Wilson loop is satu- 
rated by the quadratic (Gaussian) FC, which is proved by Casimir scaling on 
the lattice and understood as an expansion in powers of (crA)^ = 0.04 ^ 1. 

Thus all observables can be expressed via Gaussian FC with accuracy of 
few percent. On the other hand we can calculate FC as gluelump Green's 
functions, and the latter are computed again in terms of FC (or simply in 
terms of a, which is integral of FC). Thus one obtains the closed scheme 
of calculation and approximations, which is subject to the selfconsistency 
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check. This analysis reveals an interesting picture, explaining why gluonic 
condensate is relatively weak (as compared with hadron mass scale), and 
supporting quantitatively the idea of positive G2 and hence negative vacuum 
energy density, yielding np stable vacuum. Finally, the idea of QCD as 
selfconsistent theory defined by the only one mass scale, is realized in this 
FC approach and the first connection between a and Aqcd, exemplifying 
dimensional transmutation, is obtained in f lHT]) . 

The approach, discussed above, is universal and can be applied both to 
spectrum and scattering, in particular to high-energy processes, where np 
contributions are important. In all cases the problem reduces to the ground 
state or excited string dynamics, pertinent to the Wilson loop manifold with 
changeable (and integrated over) contours. Some first steps in this direction 
were done in [35l [16] . 

6 Conclusion 

The modern picture of confinement in QCD is a result of tedious efforts 
of many people during last 35 years, see |36l EZ] for reviews. The basic 
source of confinement - the nonperturbative vacuum - has passed many 
stages during this period. From constant chromomagnetic vacuum of Savvidy 
^38j . and selfdual of Leutwyler [39J, via spaghetti- vacuum of Nielsen-Olesen 
[10] . using general idea of stochastic vacuum of Ambjorn, Olesen, Peterson 
[H] and finally to the Gaussian stochastic vacuum of the Field Correlator 
Method [21 H] , discussed in this short review. Many people contributed to 
this field, suggesting models of confinement, based on different mechanisms 
(see discussion in reviews tSSllSZ]), which we could not include here, and this 
study will certainly go on. However, at present the only approach, which 
explains all details of confinement in experiment and lattice calculation, is the 
Field Correlator Method (FCM). The Gaussian correlators D and Di (and 
their colormagnetic counterparts for nonzero temperatures) are sufficient to 
describe all data up to now. The fact, that it is possible to calculate both D 
and Di selfconsistently, with the only one scale input (in addition to current 
quark masses), gives a strong hope to construct the analytic nonperturbative 
QCD ab initio, as a valuable counterpart of lattice QCD. 

The author is grateful to D.Antonov for discussions and valuable sugges- 
tions. The author acknowledges financial support of the RFFI grant 09-02- 
00629a. 
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